Imaging of Thermal Domains in ultrathin NbN films for Hot Electron Bolometers 
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We present low-temperature scanning electron microscopy (LTSEM) investigations of supercon- 
ducting microbridges made from ultrathin NbN films as used for hot electron bolometers. LTSEM 
probes the thermal structure within the microbridges under various dc current bias conditions, either 
via electron-beam-induced generation of an unstable hotspot, or via the beam-induced growth of a 
stable hotspot. Such measurements reveal inhomogeneities on a micron scale, which may be due to 
spatial variations in the NbN film or film-interface properties. Comparison with model calculations 
for the stable hotspot regime confirm the basic features of common hot spot models. 

PACS numbers: 85.25.Am, 85.25.Pb, 07.57.Kp, 74.25.Fy 



Ultrathin NbN films have recently gained much atten- 
tion for use in THz hot electron bolometer (HEB) mixers 
and fast photon counting detectors (see e.g. [H, 0, H, 3 ) ■ 
NbN HEB mixers are short microbridges (typical length 
and width of order 1 ^m and a few nm thickness) with 
contacts to antenna and dc leads. The dc current-voltage 
curve (IVC) is characteristic for microbridges with a su- 
percurrent up to a critical current Ic and a sudden jump 
to a hysteretic resistive behavior. Depending on film 
thickness, quality and geometry, typical resistances for 
this branch are around 100-1000 fi. When current bi- 
ased on the resistive branch, the bridge stays resistive 
below Ic, down to a distinct voltage from where the de- 
vice jumps back to the zero-voltage state at the return 
current 7^. When the bridge is biased at the lower sta- 
ble voltage carrying region it serves as a sensitive power 
detector for far-infrared (FIR) radiation. The basic prin- 
ciple of heterodyne mixing with such a device is the beat- 
ing of an incoming signal with an auxiliary local oscilla- 
tor (LO) signal, producing a signal at the intermediate 
frequency (IF), the difference between the frequency of 
the input signal and the LO frequency. The resulting 
IF signal lays in the 0-4 GHz range and can therefore be 
amplified with commercial semiconductor low-noise am- 
plifiers. We note that for sufficient LO power, the IVC 
becomes non-hysteretic, and the HEB mixer is typically 
voltage biased. 

NbN microbridge HEB mixers offer several very im- 
portant advantages over other mixer devices in the THz 
range. Their small size allows to pump the mixer with 
very little LO power; this is an important advantage, 
as tunable LO sources above 500 GHz have usually very 
little power output. Another advantage is due to the 
purely resistive nature of the device, which makes high- 
frequency impedance matching very easy. Because the 
detection is based on heating, the detection principle 
is not limited in frequency by specific superconducting 
properties, such as the energy gap of NbN, as in the 
case of SIS mixers. For frequencies between 1.2 and 
10 THz, HEB mixers are therefore offering currently the 
lowest noise. However, current NbN HEB mixers suffer 



from a relatively low IF bandwidth. This is due to a 
limited detection speed, which is often found to be be- 
low 2 GHz. For THz applications, IF bandwidths up to 
10 GHz are however desirable in radio astronomy. The 
detailed physics of the involved time constants, which ul- 
timately limit the IF bandwidth, is therefore a subject of 
ongoing research in many different laboratories. 

While first results for these applications are promising 
0,0], modelling of the devices is quite complex. As a con- 
sequence, the impact of fabrication and specific material 
parameters on device performance remains unclear, and 
optimization is slow and based on trial and error. Beyond 
first lumped element models , more elaborated models 
have been proposed [1, Q , which describe the physics of 
the microbridges in a spatially resolved manner. Most of 
the numerical and analytical spatially distributed models 
are based on assumptions first made for superconduct- 
ing bridges by Skocpol, Beasley and Tinkham [l3|. In 
this approach the resistive behavior of a superconduct- 
ing microbridge is modelled through the thermal equilib- 
rium between a normal conducting and therefore dissi- 
pating hot spot and the substrate. The lateral thermal 
conduction and the cooling to the substrate determine 
the detailed shape of the hot spot. The general solution 
within this analytical frame work is a symmetrical hot 
spot geometry centered between the contacts. Although 
some convincing results, such as modelling of rf-pumpcd 
IVCs and bias dependent conversion gain curves could 
be obtained for the current NbN HEBs within the exist- 
ing distributed models, the existence of a single centered 
hot spot geometry for these devices has not been directly 
shown so far. 

Over the last two decades, low-temperature scanning 
electron microscopy (LTSEM) has been used to pro- 
vide local information (on a micron scale) on various 
properties of superconducting thin films and Joseph- 
son junctions, such as the spatial distribution of the 
transition temperature Tc and critical current density 
ic , or on Josephson vortices in long junctions (Til. [T3 . 
Il3j . and on Abrikosov vortices [T^. fig and supercur- 
rent distribution [l^ in SQUID washers. Very early. 
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FIG. 1: (a) NbN HEB layout: W = 9.6 ^tm wide NbN strip, 
partially covered by Au pads with L — 4 /im separation, (b) 
IV characteristic measured at T ~ 5 K in the LTSEM. Open 
diamond and circle show bias points for LTSEM images shown 
in Fig. 2(b) and Fig. 3, respectively 



LTSEM has been applied to proof the concept of hot 
spo t formation in long and thick superconducting bridges 
[l3,[i&,19.]. In this paper we show that LTSEM can also 
be applied to investigate the thermal structures of mi- 
crobridges from ultrathin NbN films with much smaller 
geometries, approaching those of HEBs. 

The samples which we investigated were fabricated as 
follows: Thin NbN films were deposited from a 4 inch Nb 
target by 13.56 MHz (240 W) rf magnetron sputtering on 
a 2 inch fused quartz substrate at room temperature in 
a 0.852 Pa N2/Ar/CH4 atmosphere (gas flow: 2.7, 46 
and 0.6 seem, respectively). This process yields typically 
Tc = IIK for i = 5nm thick films and = 15 K for 
t > 20 nm. The NbN films were patterned by reactive ion 
etching in a CF4/O2 mixture to form long microbridges 
of width W. Subsequently. 50 nm thick Au contact pads 
were formed on top of the NbN bridges, with a separa- 
tion, which defines the length L of the active region of the 
NbN HEB, as shown in Fig.[T]^a). Typical device geome- 
tries vary from L xW = 0.5 /im x 4 /im to 5 /im x 10 /im. 
We present results from a device with t — 5 nm, L = 4 /im 
and W = 9.6 /im, which shows a clearly hysteretic IVC at 
T w 5K [c.f. Fig.mb)]. We note that all data presented 
here have been obtained with current bias. 

For imaging by LTSEM, the sample was mounted on 
a liquid He cooled stage and operated at a temperature 
T Ri 5K. The local perturbation by the focused elec- 
tron beam (e-beam) induces an increase in temperature 
5T{x,y) on the sample surface in the (x,?/) plane, on 
a length scale of approximately 1 /im, which determines 
the spatial resolution of this imaging technique, and with 
a maximum local increase in temperature AT(a;o,?/o) 
of a fraction of 1 K, centered on the beam spot posi- 
tion (a;o,yo) on the sample surface [ll|- For modelling 
the local perturbation we have used a combination of 
Monte-Carlo calculations for the beam energy deposi- 
tion and finite element modelling (FEM) of the cor- 
responding thermal plume. However, simplified semi- 
analytical models turned out to give very similar results 
[1^ . The local change in T may change global properties 
of the bridge, e.g. the voltage V across the current-biased 
bridge. This voltage change 5V depends on the e-bcam 



position (a;o,?/o) a-nd thus can be recorded to obtain a 
5V{xQ,yQ)-\ma,ge. To improve the signal-to-noise ratio, 
we use a beam-blanking unit operating at 5 kHz and de- 
tect 5V with a lock-in amplifier. 

In order to create a voltage drop along the microbridge, 
a continuous domain (hot spot) of normal conducting film 
across the bridge is required. Such a domain can be gen- 
erated either by resistive heating due to the bias current 
lb > Ir or, for lower bias currents, by the combined effect 
of e-beam and current heating. It is therefore possible to 
distinguish two different experiments: 

(A) e-beam-induced generation of (unstable) hot spot: If 
biased below Ir , the device may switch to a resistive state 
upon e-beam irradiation due to generation of an unstable 
hotspot (i.e. the hotspot disappears when the e-beam is 
turned off). The lock- in detected voltage signal corre- 
sponds to the voltage induced by the (unstable) hotspot. 

(B) e-beam induced growth of (stable) hotspot: If biased 
slightly above Ir in the resistive state, a stable hotspot is 
generated, which size may be altered upon e-beam irradi- 
ation. The lock-in detected voltage signal corresponds to 
the beam spot position-dependent small change in volt- 
age drop due to the beam-induced extension of the oth- 
erwise self-sustained hotspot. 

We first discuss imaging mode (A). Figure [2l^a) shows 
IVCs recorded during e-beam irradiation (for various 
values of the beam current lei and fixed beam voltage 
Uei = 20 kV) at a fixed position on the center of the de- 
vice, and for comparison, the IVC without irradiation. 
With increasing I^i the critical current is reduced, while 
the normal resistance remains almost unchanged. Un- 
der sufficiently strong e-beam irradiation (/e;~0.2 nA) the 
hysteresis in the IVC vanishes. The general shape of the 
IVCs with c-beam irradiation is surprisingly similar to 
IVCs of rf pumped devices This is a strong indication 
that a general mechanism, which is likely to be of ther- 
mal nature, describes the physics of the IVCs for very 
different types of energy input. 

Figure [2jb-d) shows LTSEM voltage images recorded 
at slightly sub-critical current bias {h-Ir) for fixed Ub 
and different values of lei- These images reflect the lo- 
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FIG. 2: (a) Non-hysteretic IVCs for e-beam irradiation on the 
center of the device {hi from 0.22 to 0.86 nA; Uei = 20 kV); 
unirradiated hysteretic IVC is shown for comparison (arrows 
indicate sweep direction of /;,). (b-d) LTSEM voltage images 
{lb ~ 95 /iA, i.e. just below Ir, Uei = lOkV) for different hi'- 
(b) 0.24 nA, (c) 0.26 nA, (d) 0.28 nA. The rectangles mark the 
boundaries of the sample [c.f. dashed box in Fig[lla)]. 
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FIG. 3: (a) LTSEM voltage image at = 110 fiA (stable hot 
spot condition); lei — 10 pA, Uei = 5kV. Dark areas indicate 
high voltage responsivity to e-beam irradiation, (b) Line- 
scan of LTSEM voltage signal (solid curve) along the dashed 
line in (a) and comparison with model calculation using 2- 
dimensional FEM techniques (dashed curve). Inset indicates 
the transition between superconducting (S) and normal (N) 
region at y = ±1.4 /im. 



cal sensitivity to e-beam irradiation for triggering of an 
unstable hot spot. For lei = 0.24 nA (b), no hot spot 
is induced at all. With an increase to lei ~ 0.26 nA (c) 
a small spot on the image appears: only when the e- 
beam is applied within this spot a voltage is induced. As 
expected from thermal and electrical symmetry consid- 
erations, the spot is centered with respect to the x-axis. 
It is also centered with respect to the y-axis (along the 
bridge), which can be easily explained by the additional 
cooling effect provided by the gold pads. With further 
increase of lei to 0.28 nA, the spot increases rapidly [see 
Fig. [2l^d)] . An analysis of the shape of the spots leads to 
the estimate of the spatial resolution of ~ 0.4 fj,m for this 
type of measurement. 

For imaging mode (B), (/;, > J^), i.e. in the area of 
a stable bias current driven hot spot, the voltage across 
the sample increases slightly under irradiation due to the 
following effect: The beam-induced increase in tempera- 
ture causes the hot spot to grow; thus an increased re- 
gion of the NbN film becomes resistive, and the normal 
resistance R is increased. This mechanism is most effec- 
tive for an e-beam position close to the superconduct- 
ing/normal conducting (S/N) boundary. The result is a 
characteristic double ridge picture, with the maximum of 
the signal indicating the location of these lateral bound- 



aries, as shown in Fig.[3lja). However, the LTSEM signal 
also reveals a clear variation along the S/N boundaries. 
This may be due to inhomogeneities in the film quality 
such as thickness or composition, or due to local vari- 
ations in the thermal coupling to the substrate or con- 
tact pads. The detected inhomogeneities will result in 
an asymmetric current density and may adversely affect 
device performance. 

We modelled the device under the stable hotspot bias 
conditions as in imaging mode (B) by a 2-dimensional 
finite element calculation, taking into account the heat 
transfer to the substrate as balanced by the heat pro- 
duction by the bias current dissipation. The e-beam in- 
duced change of the substrate surface temperature was 
modelled using Monte-Carlo techniques to compute the 
energy input and a finite element model in the symme- 
try plane of the e-beam axis to derive the temperature 
resulting from this energy input. For the parameters 
(/{,. lei, Uei) and device geometry used in the exper- 
iment [c.f. FiglS], the model predicts a superconduct- 
ing (S)/normal conducting (N) transition at a position 
of y = ±1.4 ^m, as indicated in Fig[3l^b). The derived 
voltage signal SV{y) for a linescan along the y-direction 
[c.L Fig.[3jb)] reproduces nicely the characteristic double 
peaked shape as observed experimentally, and matches 
also very well the general form of the measured curve. 

In conclusion, we demonstrated that LTSEM is a use- 
ful tool to investigate thermal domain (hotspot) forma- 
tion in ultrathin superconducting films. Our experiments 
show that for current bias below Ir, e-beam irradiation 
induces a hotspot, depending on deposited beam energy, 
beam spot position and sample inhomogeneity. For bi- 
asing above 7^, i. e. in the regime of a stable hotspot, 
comparison of measurements with 2-dimensional mod- 
elling shows that the classical hotspot model does apply. 
Our investigations also show that inhomogeneities on a 
micron scale can be detected. Further investigations are 
required to clarify the nature of such inhomogeneities and 
their impact on device performance. Furthermore, it will 
be interesting to include RF pumping to investigate the 
hotspot formation under typical working conditions of 
superconducting hot electron bolometers. 

D.D. gratefully acknowledges support from the Evan- 
gelisches Studienwerk e.V. Villigst. This work was sup- 
ported by the EU FP6 Program RADIONET (AM- 
STAR). 



[1] Y. P. Gousev et al, J. Appl. Phys. 75, 3695 (1994). 
[2] J. Kawamura et al., J. Appl. Phys. 80, 4232 (1996). 
[3] T. Lehnert et al, J. Appl. Phys. 83, 3892 (1998). 
[4] G. N. Gol'tsman et al, Appl. Phys. Lett. 79, 705 (2001). 
[5] S. Cherednichenko et al, Physica C 372-376, 427 
(2002). 

[6] M. C. Wiedner et al, A&iA 454, L33 (2006). 
[7] B. S. Karasik, A. I. Elantiev, Appl. Phys. Lett. 68, 853 
(1996). 

[8] D. W. Floet et al, Appl. Phys. Lett. 74, 433 (1999). 
[9] H. F. Merkel et al, IEEE Trans. MTT 48, 690 (2000). 
[10] W. J. Skocpol et al, J. Appl. Phys. 45, 4054 (1974). 



[11] J. R. Clem, R. P. Huebener, J. Appl. Phys. 51, 2764 
(1980). 

[12] R. P. Huebener, Rep. Prog. Phys. 47, 175 (1984). 
[13] R. Gross, D. Koelle, Rep. Prog. Phys. 57, 651 (1994). 
[14] D. Doenitz et al, Appl. Phys. Lett. 85, 5938 (2004). 
[15] R. Straub et al, Appl. Phys. Lett. 78, 3645 (2001). 
[16] D. Doenitz et al., Phys. Rev. B 73, 064508 (2006). 
[17] R. Eichele et al, Appl. Phys. Lett. 38, 383 (1981). 
[18] R. Eichele et al, Z. Phys. B 48, 89 (1982). 
[19] R. Eichele et al, J. Low. Temp. Phys. 52, 449 (1983). 



